Field-Dependent Critical Current in Type-II Superconducting Strips: Combined 
Effect of Bulk Pinning and Geometrical Edge Barrier 
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Recent theoretical and experimental research on low-bulk-pinning superconducting strips has re- 
vealed striking dome-like magnetic-field distributions due to geometrical edge barriers. The observed 
magnetic-flux profiles differ strongly from those in strips in which bulk pinning is dominant. In this 
paper we theoretically describe the current and field distributions of a superconducting strip under 
the combined influence of both a geometrical edge barrier and bulk pinning at the strip's critical 
current I c , where a longitudinal voltage first appears. We calculate I c and find its dependence upon 
a perpendicular applied magnetic field H a - The behavior is governed by a parameter p, defined 
as the ratio of the bulk-pinning critical current I p to the geometrical-barrier critical current I 3 o- 
We find that when p > 2/n and I p is field-independent, I c vs H a exhibits a plateau for small H a , 
followed by the dependence I c — Ip oc H^ 1 in higher magnetic fields. 



The combination of a geometrical edge barrier and 
bulk pinning recently has been shown to strongly affect 
the properties of low-dimensional superconductors (thin 
films, single crystals, and tapes with high demagnetiz- 
ing factors) placed in either a perpendicular magnetic 
field 1-5 or a transport-current-carrying state 6-9 . While 
most experimental studies of the field dependence of the 
critical current I c are being interpreted solely on the ba- 
sis of bulk-pinning theory (see for example 10- ), a num- 
ber of works 6,8 ' 15,16 have shown that a geometrical edge 
barrier (or surface barrier) may strongly affect I c . In this 
paper we study the combined effect of a geometrical edge 
barrier and bulk pinning upon the magnetic field depen- 
dence of I c for type-II superconducting strips. We shall 
show how the dependence of I c upon H a is controlled by 
the parameter p = I p /I s o, where I p is the bulk-pinning 
critical current in the absence of a geometrical edge bar- 
rier, and Iso is the geometrical-barrier critical current in 
the absence of bulk pinning. 

We consider a superconducting strip of thickness d 
(\y\ < d/2) and width 2W (\x\ < W) centered on the 
z axis. We assume that d is less than the London pen- 
etration depth A and that W is much larger than the 
two-dimensional screening length A = 2X 2 /d. The strip 
is subjected to a perpendicular applied magnetic field 
H a = (0, H a ,0), and it carries a total current I in the 
z direction described by a spatially dependent sheet cur- 
rent density K.(x) = 3d = [0, 0, K z {x)]. We wish to deter- 
mine the current-density and magnetic-field distributions 
at the critical current at which a steady-state flux-flow 
voltage appears along the length of the strip. For a strip 
containing no magnetic flux, K z (x) is the sum of two 
contributions, 8 
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the Meissner-state current density induced by the applied 
field H a . The divergences in Eqs. (1) and (2) at |a;| = W 
are cut off when x is within A of the edge. 

To account for the edge barrier, we assume that vor- 
tices nucleate and enter the superconductor when K z 
at either sample edge reaches the value K s — j s d at 
which the barrier is overcome. For an ideal edge, j s is 
equal to the Ginzburg-Landau depairing current density 
jgz, 18,19 , but for an extremely defected edge, j s may be- 
come negligibly small. When H a = 0, the sheet cur- 
rent at both edges is approximately K] z w I /ir\/2WA, 
such that the edge-barrier critical current in zero exter- 
nal magnetic field is I s q « ttK s \/2WA. When H a > 0, 
the net sheet current at x = W is approximately K z sa 
(I + 2nH a W)/(TrV2WA), and the edge-barrier critical 
current becomes I s (H a )/ I s q = 1 — h for small h, where 
h = H a j (I s0 /2ttW). This result is essentially the same as 
that found in Ref. 8 for the critical current in low applied 
magnetic fields for bulk-pinning- free strips. 

We next account for bulk pinning, characterized via 
a bulk-pinning critical sheet current density, K p = j p d, 
such that the critical current in the absence of an edge 



barrier is /„ 



2K p W. 



We first consider the case 
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of relatively weak bulk pinning when I p < (2/tt)I s q, 
i.e., p < 2/ir. In low fields (0 < H a < Hd, region I 
of Fig. 1), vortices nucleate on the right-hand side at 
x = W when / slightly exceeds I S (H). As long as 
K z (x) — Ki z (x) + K az (x) exceeds K p , these vortices 
are driven entirely across the strip, traveling rapidly 
(speed v governed solely by the force-balance equation 
[K z (x) — K p ](f>o — rjvd, where r\ is the viscous drag coef- 
ficient), and annihilating with their images on the op- 
posite side of the strip. The critical current is then 
I c {H a ,p) = I s (H a ), and the normalized critical current 



the Meissner-state current density generated by the ap- 
plied current /, and 



i c {h,p) = I c {H a ,p)/I s0 = 1 - h. 



(3) 
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However, K m i ni the minimum value of K z {x), decreases 
with increasing H a and reaches K p at / = I s (H a ) when 
H a = H d , where 



h d = H d /(I s0 /2nW) = hi - (ttp/2) 2 } 



(4) 




FIG. 1. Behavior at the critical current vs reduced field 
h and bulk pinning parameter p. In region I, the strip is 
vortex-free [H y (x,0) = 0] and the sheet current density is ev- 
erywhere above the bulk-pinning critical value [K z (x) > K p ]. 
In II, there is a vortex-free zone where K z (x) > K p on the 
right side of the strip and a vortex dome (dot-dashed curve in 
inset) where K z (x) — K p (solid curve in inset). In III, there 
are three zones: two vortex-free zones on either side of a vor- 
tex dome. In IV, there are four zones: a vortex dome where 
H y (x,0) > 0, an antivortex dome where H y (x,Q) < 0, and 
two vortex-free zones where K z (x) > K p near the edges. The 
curve hd(p) (solid) separates regions I and III, hi(p) (dashed) 
separates III and IV, and h,2 (p) (dotted) separates II and III. 

When p < 2/tt and h > h d or when p > 2/tt, i.e., for h 
and p outside region I of Fig. 1, nucleated vortices stop 
inside the strip, and dome-like flux distributions occur 
at the critical current. To determine the critical current 
with domes present, one must first calculate the vortex 
(and antivortex) density n(x) = fioHylx, 0)/<po inside the 
dome, where K z (x) = K p , and the sheet current density 
K z {x) outside, where n{x) = 0. We have obtained these 
mathematically by using the Cauchy integral inversion 
method 2 ' 4 ' 7,9 ' 20 to invert the Biot-Savart law. However, 
we shall use the method of complex fields 21 to give a 
physical interpretation of the mathematical results. 

Following Ref. 21 , we express the two-dimensional field 
distribution as an analytic function H(() — H y (x,y) + 
iH x (x, y) of the complex variable Q — x + iy, such that 
the Biot-Savart law becomes 



H(Q = H a + 
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The inversion procedure yields H(() of the following form 
at the critical current I c (H a ): 

H \0~ ( C 2 _ W 2y/2 [H a + —Q(a,b,Q\, (6) 



where 
Q(a,b,Q 



= fb y/W^f du 

Ja (C-u)^{u-a)(b-u) 

2{W+a) r jW-Q jjf iC+W)^-^ 

L (C-a) ±1 \(C-a)(W+b)^J 



(7a) 



y/(W-a)(W+b) L TC^T i±V (C-a)(W+b) . 

+ n(^,g)], (7b) 



and 



2 _ 2W(b-a) 



[W -a)(W + b)' 
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In Eq. (6), the term proportional to H a is simply the com- 
plex field describing the Meissner-state response to the 
applied field H a of two parallel superconducting strips 8 ' 22 
(-W < x < a and b < x < W). The term propor- 
tional to K p is the complex field describing the image- 
current response 23 of the two strips to currents K p du 
summed over the region a < u < b. We have evalu- 
ated the integral in Eq. (7a) in terms of complete ellip- 
tic integrals of the third kind 24-28 H(n,k), where n is 
called either the characteristic or the parameter, and k is 
called the modulus. Equation (7b) can be used to evalu- 
ate H y (x, 0) = RcH(x) and K z (x) = -2ImH(x + ie). 

For p > 2/n and small values of h, i.e., for h and p in 
region IV of Fig. 1, the vortex distribution at the critical 
current can be described as a double dome, consisting 
of a vortex dome adjacent to an antivortex dome (see 
inset). Just above the critical current, vortices nucle- 
ate at x = W(x' = 1), where K Z (W — A) = K s , move 
rapidly to the left through an otherwise vortex-free region 
(b < x < W), and then move slowly to the left through a 
vortex-filled region (the vortex dome), a + < x < b. An- 
tivortices nucleate at x = —W, where K z {— W+K) = K s , 
move rapidly to the right through an otherwise vortex- 
free region (— W < x < a), and then move slowly to the 
right through an antivortex-filled region (the antivortex 
dome), a < x < 6_. Vortices and antivortices annihilate 
where the two domes meet at x = b- = a + . 

Two equations must be solved simultaneously for a 
and b at the critical current I c (H a ) for known values of 
h and p in region IV of Fig. 1. One condition is that 
K Z {W - A) = K s , which yields from Eqs. (6) and (7b) 



y/(l-a')(l-V)h- 



(1 + 
= 1 



b' tt / b' —a 
1+b' 1+6' 



q)p 
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where we use the normalized quantities a' = a/W and 
b' = b/W. The other condition, that K z (-W + A) = K s , 
yields 
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- y/(l + a>)(l + b>)h + (1 - b')^^H( b T E^,q)p 

= i. ' (io) 

Expansion of Eq. (5) for large £ yields -ff(C) = #a + 
I/2tt( + 0(1/( 2 ). Expanding Eqs. (6) and (7a) in powers 
of ( and making use of Eq. (9), we obtain the normalized 
critical current, i c = I c /I s o- 

Uh,p) = 

CV l! 2V(l-a')(l-6') 

+ 1 V^K 1 - a ')E(g) + (1 + a')K(«)]p, (11) 

where a' and 6' are determined from Eqs. (9) and (10) 
for the desired values of h and p. 

The double-dome vortex-antivortex distribution (re- 
gion IV of Fig. 1) occurs at i c only for h in the range 
< h < hi for p > 2/tt. Here, h\(p) is the low- 
est value of h that makes H y (x,0) > in the region 
a < x < b. Thus, one of the equations determining h\ is 
H a + (K p /2n)Q(a,b,a + e) = [see Eq. (6)], which yields 



(b'-a')y/(l-a')(l+b') 

[(l-a')(l + b')E(q)-(l + a')(l-b')K(q) 

-(b'-a')(l-b')Tl(^, q )]p = 0. (12) 

hi can be determined for a given value of p as the value of 
h when Eqs. (9), (10), and (12) are simultaneously solved 
for a', b', and h. 

For known values of h and p in region III of Fig. 1, the 
left and right boundaries of the vortex dome a' and b' are 
determined by simultaneously solving Eqs. (9) and (12); 
Eq. (12) also gives the condition that dK z {x)/dx = at 
x = a. Once a 1 and b' are found, Eq. (11) again can be 
used to calculate the critical current. Just above the crit- 
ical current, vortices nucleate at x = W , move rapidly to 
the left through the vortex-free region b < x < W, move 
slowly to the left through the vortex dome a < x < b, 
escape from the dome, and finally move rapidly to the 
left through the vortex- free region — W < x < a. 

For increasing values of h, the left boundary of the 
vortex-hlled region moves closer to the left edge of the 
strip; a becomes equal to — W + A when h = h%, which 
can be determined for a given value of p as the value 
of h when Eqs. (9) and (12) are numerically solved for 
b' and h, taking a' = -1 + A/W. (For Figs. 1 and 2, 
A/W = 0.01 was assumed.) For h > h2, the field and 
current distributions and i c can be calculated with good 
accuracy by simply setting a = W in Eq. (6). The com- 
plex field in region II of Fig. 1 is then 

H M = §^- 2 [H a + ^Q(b,0], d3) 

where 

W0 = J!^^_ (14a) 

= 2sinh- V ^-2sinh-7{f±8g^. (14b) 



The condition K Z (W — A) = K s , which determines 
b' = b/W at I c for h and p in region II of Fig. 1, be- 
comes 

^(l-b'^h+psmh- 1 y/(l + 60/(1 -V)] = 1 (15) 

instead of Eq. (9), and the normalized critical current 
can be expressed as 



t c (h,p) = -^2(1 -b') + 1^2(1 + 60 (16) 
instead of Eq. (11). 
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FIG. 2. i c (h,p) (critical current normalized to I s o) vs re- 
duced field h for fixed values of the bulk pinning parameter p. 
The solid straight line and solid circles denote values of i c in 
region I, where h = hd{p); the dashed curve and solid squares 
show i c at h — hi(p); and the dotted curve and solid triangles 
show i c at h — fe(p). For p < 2/tt, i c decreases linearly with 
h [Eq. (3)] up to hd and then decreases more slowly in regions 
III [Eq. (11)] and II [Eq. (16)]. The bold curve shows i c for 
the special case of p = 2/tt. For p > 2/tx, i c increases by a 
few percent in the double-dome region IV and then decreases 
more gradually in regions III and II. In all cases, i c asymptot- 
ically approaches p for large h (short dotted lines along the 
right side of the figure) . 

The reduced critical current i c as a function of h, calcu- 
lated from Eq. (11) or (16), is shown in Fig. 2 for various 
values of p. For p = and h > 1/2, we find i c = l/4/i, as 
obtained in 6,8 for pin-free strips with an edge barrier. In 
the opposite limit, p 3> 1, we obtain i c w p, as expected 
for bulk-pinning-dominatcd behavior. For fe> p, we see 
from Eq. (15) that b' approaches 1, and Eq. (16) yields 
i c w p + l/4/i. A generic feature of Fig. 2 is the plateau 
in i c vs h to the left of the dashed curve in region IV; ac- 
tually, i c increases by a few percent as h increases from 
to hi . This increase is due to a significant decrease in the 
width b — a over which K z (x) is restricted to K p . This 
effect is partially compensated by a change in shape of 
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K z (x) in the vortex-free zones [e.g., dK z (x)/dx — at 
x = a at h = hi]. Field-dependent critical current densi- 
ties j c (H a ) have been found experimentally in Refs. 10-14 , 
but the behavior was interpreted solely in terms of bulk 
pinning. 

In agreement with earlier work, 9 ' 29 our results show 
that the critical current of a strip is not a simple su- 
perposition of currents / s and I p , as was suggested in 
Rcfs. 16 ' 30 ' 31 . Only in the limit h > p is it possible to 
express the critical current as I c {H a ) = I p + I s (H a ). 

We have assumed here that K p is a constant. Because 
of the nonlocal current-field relation [Eq. (5)], it would 
be a challenging task to find I c (H a ) when K p depends 
upon the local magnetic field H y (x, 0). 

In summary, we have solved for the field dependence 
of the critical current density in a superconducting strip 
accounting for both bulk pinning and a geometrical edge 
barrier, and we have developed a procedure for finding 
the magnetic-field and current-density distributions in- 
side the strip at the critical current. In the presence of a 
strong edge barrier, we have found strong field dependen- 
cies of the critical current. Such effects should be taken 
into account when interpreting experimental critical cur- 
rents in low and moderate magnetic fields. 
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